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Abstract 

The Boltzmann type kinetic equation for solitons in Nonlinear Schrodinger equation has 
been constructed on the base of analysis of two soliton collision. Possible applications 
for Langmuir solitons in plasma and solitons in optic fibers are discussed. 
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1. Nonlinear Schrodinger equation (NSE) is one of the most popular equations for 
describing nonlinear phenomena in condensed matter and plasma. Such effects as small 
amplitude localized waves in magnets, localized excitations in quasi-one-dimensional 
biologic systems and optic fibers are described by NSE (see for example [1-3]). In spite 
of wide field of applications, kinetic properties of solitons in NSE as well as the another 
models were investigated much less then kinematic and dynamic properties of solitons. 
Some aspects of kinetic behavior of solitons in Korteveg-de-Vries, 0^ and Sine-Gordon 
equations have been investigated in [4-8]. Fokker-Plank equation for the distribution 
function of soliton in NSE damping and fluctuations were considered in [9]. In present 
paper the kinetic equation for solitons in NSE following the approach suggested in [10] 
is proposed. Possible applications for solitons in plasma and optic fibers are discussed. 

2. NSE in dimensionless form for slow varying nonlinear envelope is written as: 

Construction of kinetic equation for solitons needs analysis of the process of soliton- 
soliton interaction. It is well known that solitons in exactly integrable models do not 
change there forms and velocities and additional shifts of there coordinates Ax and 
phases A0 appear as the result of their interaction [1-3]. For NSE the shifts Ax and Acf) 
can be written as [1] : 

= Cr'ln|Zi,2| A0i = 2ar^Zi,2, i = l,2 (2) 

^ ^ 2^(1 + (2) - (ri- V2) 

Here Q - is dimensionless parameter characterizing width and amplitude of soliton, Vi 
-is soliton velocity. 

To simplify the problem we will analyse the case of low density of soliton, therefore 
we will consider the case when shift of soliton coordinate much more then soliton width. 
In the region of parameters: 

4(Ci - < (^1 - v^r « 4(Ci + C2)^ (3) 

one can obtain that 

Axi > Cr\ > C2~' (4) 

and 

\A(f)i\ ^ \A(f)2\ ^ 1, (5) 

that gives us the possibility to follow only for the evolution of soliton positions. 

Let us introduce the distribution function of solitons f{x,p, t). To construct collision 
integral we follow [10] and apply Boltzmann approach taking into account that in con- 
trast to usual particles solitons change there positions. Also it is necessary to consider 
that the interaction between solitons is absent when the distance become much more 
then soliton width. Let us analyse the number of solitons arriving to the point {x,pi) 
and leaving this point due to two soliton collisions for example in the case P2 > Pi > 0. 
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Due to following process: 



{x,pi) {x + Axi,pi) 

(X,P2) ^ (X- AX2,P2) (6) 



the number of solitons leaving the point x, pi can be written as: 

Ar(-) ^dx [ \vi -V2\fix,pi,t)f{x,p2,t)dp2. (7) 

The number of solitons arriving at the point {x,pi) and corresponding to following 
process 

{x - Axi,pi) {x,pi) 

{X- Axi,P2) ^ {X- Axi- AX2,P2) (8) 

can be presented as: 

N^+^^dx r \vi-V2\f{x-Axi,pi,t)f{x-Axi,p2,t)dp2. (9) 

Jv-i 

From formulacs (7) and (9) one can write soliton-soliton collision integral in the case 
P2 > Pi > 0. Considering all possible ratio for pi and p2 soliton-soliton collision integral 

can be written as: 

/oo 
l^^i - 'V2\{f{x - Axi,pi,t)f(x - Axi,p2,t) - f(x,pi,t)f(x,p2,t)}dp2. (10) 
-oo 

In the approximation 

\Axi\ pa |Aa;2| ~ \Ax\, Ax = Q^^signivi — V2) ln( Y (11) 

Vi - V2 

collision integral coincide with one obtained for Sine-Gordon equation in [10]. Assuming 
that f{x,p,t) is slowly varying in scales comparable to Ax and expanding f{x,p,t) in 
powers of Ax and keeping the leading terms we can rewrite the expression for C in the 
form: 

d 1 d"^ 

where 

u{x,pi,t)^ \vi -V2\Axf{x,p2,t)dp2, (13) 
J —00 

describes the renormalization of soliton velocity and 

/oo 
\Vi-V2\{Axff{x,P2,t)dp2, (14) 
-00 

describes the solitons diffusion process. 
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The kinetic equation can be written as: 

df{x,p\^t^ d 1 d'^ 

dt ^ = l^-Q^T^{^^Pi^'t)f{^^Pu't). (15) 

From the kinetic equation (15) it is easy to show, that soliton-soliton collisions lead to 
the entropy production in the soliton gas in the case of nonuniform distribution function 
of soliton in coordinate space (see [10]). Applying standard methods from (10) it is easy 
to derive transport equations and calculate the relaxation time (see [10,11]). 

3. Let us consider some apphcations that result from kinetic equation (15). NSE 
appears in analysis of nonhnear electron plasma waves [12-15]. The conditions of solitons 
creation [16] and collapse [17,18] and interpretation of Langmuir turbulence in soliton 
terms have been analysed in the frame of NSE (see [13-15]). 

On the base of proposed approach it is possible to estimate the relaxation time Tgsfor 
soliton gas. Indeed this estimation can be written as: 

— {qxofnvT, (16) 

where xq - is soliton size, q is characteristic scale of inhomogeneity in soliton gas, n- 



is soliton density, vt = yT/m - is thermal velocity of solitons, T -is temperature of 
solitons and m is its mass. 
Here 

where Vg - is group velocity of the high frequency wave, ou and k - its frequency and wave 
vector, ijj -is slow varying nonlinear envelope. 

To come nearer from proposed scheme to real situation it is necessary to take into 
account two facts: the interactions which destroy the integr ability of the system are 
always existed and one dimensional soliton is unstable under the influence of 2D and 3D 
perturbations. In the nonintegrable system solitons collide with momentum changing, 
but in the case close to completely integrable model it is not difficult to estimate the 
relaxation time 1/t*^ ~ AE/E, where AE/E is relative energy change due to coUision. 
Therefore we have two steps of the relaxation: first one deal with shifts of the soliton 
positions and second deal with momentum exchange [10,19]. The characteristic time r„„ 
of soliton instability in 2D space is finite and r^^ ~ (see [14]), where L is characteristic 
size in perpendicular direction. Comparing these estimations with formula (16) it is 
possible to conclude that above considered kinetic behaviour of sohtons deal with shift of 
its positions can be realized as a intermediate regime before solitons collapse and can be 
interpreted as soliton turbulence. Soliton turbulence phenomena have been considered in 
[15] were kinetic equation for many particle distribution function of soliton like waves has 
been proposed. In fact in [15] processes with momentum exchange only were considered. 

4. Another application of NSE deal with soliton propagation in optic fibers [20,21]. 
For single mode case NSE can be written as: 
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Here X is amplitude envelope of the pulse, ^ = z / L,^ - is the coordinate along the fiber, 
Li^ - is characteristic disperse length, r = {z — vot)/vQTQ, where vq -is the group velocity. 
To is the initial duration of the pulse, the value a is definited by dispersion of group 
velocity and refraction coefficient. 

From applied point of view the problem of relaxation of solitons interacting with 
defects is one of the most actual among kinetic effects. The possible results of soliton - 
defect interaction are transmission, reflection or capturing of soliton by defect (see [22- 
25]). Dynamic properties of solitons in optic fiber interacting with dispersion-spectrum 
inhomogeneities have been studied in [26]. Obviously that the most natural way to 
analyze relaxation of solitons is to formulate soliton-defect collision integral. In general 
form it can be written as: 

/oo 
W{puP2){f{x,Put)- f{x,p2,t)}dp2. (19) 
-oo 

Concrete expression for W{pi^p2) in Sine-Gordon model was calculate in [11] for the 
case of elastic scattering of soliton by impurity. Here we estimate the relaxation time Tsd 
in the case of low density of solitons comparing with defects concentration Cj. Really 
from dimensionally consideration: 

- ~ a^. (20) 

where ^Esd is relative energy loss due to interaction with defect. 

The defects appearance can deal with fiber irradiation. As a result the characteristics 
of light and sound propagation have changed under the influence of the ionizing radiation. 
These effects have been apphed for measuring of the radiation dose [27,28]. Obviously 
that the relaxation time and kinetic coefficients of solitons depend from the intensity of 
fiber irradiation. 

Let us consider the expression for the concentration of defects in the case of 7-rays. 
The 7-rays can be generated in the processes of electron beam braking on heavy elements 
targets. The breaking radiation has continuous spectrum with the maximum energy of 
photons equal to the energy of electrons (see for example [29]). In this case [30]: 

Ci = UeK, (21) 

where t - is the time of irradiation, Je -is the current of electrons, coefficient i^=2.1 for 
Si [30] and energy of 7-rays >20Mev. 

Therefore: ^ 

"^^ " Tsd t {AEsd/E) • ^^^^ 

Consequently due to stability of soliton and low damping in optical fiber this effect 
can be applied for measurement of irradiation dose. 

This work can not been done without support of International Atomic Energy Agency, 
contract No 7996. 
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